arXiv:1503.09148vl [hep-ph] 31 Mar 2015 


Light Quarks in the Screened Dyon-Anti-Dyon Coulomb Liquid Model II 


Yizhuang Liuj^j Edward ShuryakJT] and Ismail ZahecQ 
Department of Physics and Astronomy, Stony Brook University, Stony Brook, New York 11794-3800, USA 

(Dated: April 1, 2015) 

We discuss an extension of the dyon-anti-dyon liquid model that includes light quarks in the dense 
center symmetric Coulomb phase. In this work, like in our previous one, we use the simplest color 
SU(2) group. We start with a single fermion flavor Nf = 1 and explicitly map the theory onto a 
3-dimensional quantum effective theory with a fermion that is only Uy(l) symmetric. We use it to 
show that the dense center symmetric plasma develops, in the mean field approximation, a nonzero 
chiral condensate, although the ensuing Goldstone mode is massive due to the Ua( 1) axial-anomaly. 

We estimate the chiral condensate and a, p meson masses for Nf = 1. We then extend our analysis 
to several flavors Nf > 1 and colors N c > 2 and show that center symmetry and spontaneous chiral 
symmetry breaking disappear simultaneously when x = Nf/N c > 2 in the dense plasma phase. 

A reorganization of the dense plasma phase into a gas of dyon-antidyon molecules restores chiral 
symmetry, but may preserve center symmetry in the linearized approximation. We estimate the 
corresponding critical temperature. 

PACS numbers: 11.15.Kc, 11.30.Rd, 12.38.Lg 


I. INTRODUCTION 

This work is a continuation of our earlier study [T] of 
the gauge topology in the confining phase of a theory with 
the simplest gauge group 517(2). We suggested that the 
confining phase below the transition temperature is an 
“instanton dyon” (and anti-dyon) plasma which is dense 
enough to generate strong screening. The dense plasma 
is amenable to standard mean field methods. 

While an extensive introduction to the subject can be 
found in [I], here we only mention few important points. 
The treatment of the gauge topology near and below T c is 
based on the discovery of KvBLL instantons threaded by 
finite holonomies [5] and their splitting into the so called 
instanton-dyons (anti-dyons), also known as instanton- 
monopoles or instanton-quarks. Diakonov and Petrov 
[3j suggested that the back reaction of the dyons on the 
holonomy potential at low temperature may be at the ori¬ 
gin of the disorder-order transition of the Polyakov line. 
A very simple model of a de-confinement transition has 
been proposed by Shuryak and Sulejmanpasic [3] through 
the use of dyon-antidyon “repulsive cores”. 

The dyon-anti-dyon liquid model proposed by Di¬ 
akonov and Petrov j3] was based on (parts of) the one- 
loop determinant providing the metric of the moduli 
spaces in BPS-protected sectors, purely selfdual or anti¬ 
selfdual. The dyon-antidyon interaction is not BPS pro¬ 
tected and appears at the leading - classical - level, re¬ 
lated with the so called streamline configurations, the so¬ 
lutions of the “gradient flow” equation. These solutions 
have been recently derived by Larsen and Shuryak [5]. 
Their inclusion in our work [Tj reveals a very strong cou- 
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pling of the dyons to the anti-dyons, which can however 
be effectively reduced by screening, provided the dyon 
ensemble is dense enough. 

Before turning to the main subject of this work which 
is focused on the effects of light quarks on the gauge 
topology and chiral symmetry, we will briefly mention 
some important studies for the development of our work. 
The original discovery of the KvBLL instantons [2j with 
non-trivial holonomies is the key starting point for assess¬ 
ing the role of center symmetry on the gauge topological 
structures. The second important development is the 
assessment of the quantum weight around the KvBLL 
instantons in terms of the coordinates of the instanton- 
dyons developed by Diakonov and collaborators 0 0. 
The dissociation of instantons into fractional constituents 
is similar to the Berezinsky-Kosterlitz-Thouless (BKT) 
transition in 2-dimensional CPN models [7], as has been 
advocated by Zhitnitsky and collaborators [S] , although 
substantially different in the details. 

A center-symmetric (confining) phase can be compati¬ 
ble with an exponentially dilute regime that is controlled 
semi-classically, as shown by Unsal and Yaffe [S] using 
a double-trace deformation of Yang-Mills action at large 
N on S 1 x R 3 . A similar trace deformation was used 
originally in the context of two-dimensional (confining) 
QED with unequal charges on S 1 x R ms to analyze 
center symmetry and its spontaneous breaking. This con¬ 
struction was extended to QCD with adjoint fermions by 
Unsal m, and by Unsal and others [12] to a class of 
deformed supersymmetric theories with soft supersym¬ 
metry breaking. While the setting includes a compact- 
ification on a small circle, with weak coupling and an 
exponentially small density of dyons, the minimum at 
the confining holonomy value is induced by the repulsive 
interaction in the dyon-anti-dyon pairs (called bions by 
the authors). A key role of supersymmetry is the can¬ 
cellation of the perturbative Gross-Pisarski-Yaffe-Weiss 
(GPYW) holonomy potential [15] . While this allows to 
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study deconfinement transition in the very dilute regime, 
the major subject to be studied in this work - sponta¬ 
neous chiral symmetry breaking - would still be absent, 
as its development would require an ensemble which is 
sufficiently dense. 

Let us now turn to the effects of light fermions. Key to 
these effects are topological index theorems, which relate 
the topological charge of the solitons to the number of 
its fermionic zero modes. When the ensemble of topolog¬ 
ical solitons is dense enough, the fermionic zero modes 
can collectivize and produce the so called Zero Mode 
Zone (ZMZ) which breaks spontaneously chiral symme¬ 
try. For ensemble of instantons this phenomenon has 
been studied in great detail in 1980’s and 1990’s, for a 
review see |14j . Thanks to topology, the fermionic zero 
modes are remarkably stable agains any smooth defor¬ 
mations of these objects, resisting tremendous amount of 
perturbative noise. As has been derived in the “instan- 
ton liquid model” context and many times observed in 
lattice numerical simulations, the ZMZ states with Dirac 
eigenvalues in the range |A| < 20 MeV , are crucial for the 
generation of the hadronic masses and properties, while 
being only a tiny subset of all fermionic states (typically 
of the order of 10” 4 in current lattice simulations). 

Since the instanton-dyons carry topological charge, 
they should have zero modes as well. On the other hand, 
for an arbitrary number of colors N c those topological 
charges are fractional 1 / N c , while the number of zero 
modes must be integers. Therefore only some instanton- 
dyons may have zero modes. For general N c and general 
periodicity angle of the fermions the answer is known 
but a bit involved. For SU(2) colors and physically anti- 
periodic fermions the twisted L dyons have zero modes, 
while the usual M-dyons do not. 

Recent investigations by Shuryak and others nani 
have shown that light fermions disorder in ensembles 
composed of interacting dyons and anti-dyons. Their nu¬ 
merical analyses support the accumulation of zero virtu¬ 
ality quark states leading to the spontaneous breaking of 
chiral symmetry. The effects of light quarks on deformed 
supersymmetric models where recently addressed in Em¬ 
in this work we would like to follow up on our study 
in |T] by allowing for light quarks in the dense center 
symmetric phase of the dyon-anti-dyon Coulomb plasma. 
The word “dense” is key here, as it justifies the use of a 
mean-held analysis in the characterizing the spontaneous 
breaking of chiral symmetry and the formation of a chiral 
condensate. As our interest is now in the light quark dy¬ 
namics, we will only enforce the strong Coulomb correc¬ 
tions at the constraint level. One of the chief achievement 
of this work is to show how the induced chiral effective 
Lagrangian knows about confinement. In particular, we 
detail the interplay between the spontaneous breaking of 
chiral symmetry and center symmetry. 

In section 2 we detail the color SU(2) version of the 
model for Nf = 1. By using a series of fermionization and 
bosonization techniques we show how the 3-dimensional 
effective action for the liquid can be constructed to ac¬ 


commodate for the light quarks. In section 3, we show 
that the ground state solution supports both center sym¬ 
metry and chiral condensation. In section 4 we detail the 
flavor spectrum in terms of the sigma meson, the eta' 
meson which is shown to be anomalous. In section 5, 
we explore the effects of molecular pairing of dyons and 
anti-dyons induced by the light quarks near the tran¬ 
sition temperature and their effect on the formation of 
the chiral condensate and center symmetry. In section 
6 we briefly extend the model to include many colors 
and flavors and show that in the dyon-anti-dyon liquid 
with light quarks, the restoration of chiral symmetry oc¬ 
curs simultaneously with the loss of center symmetry for 
x = Nf /N c > 2. An estimate of the transition tempera¬ 
ture from the center symmetric to non-symmetric phase 
is made. Our conclusions are in section 7. 


II. EFFECTIVE ACTION WITH FERMIONS 
A. General setting 

Since this is the second paper of the series, we will 
keep all notations consistent with the first paper [T[ which 
should be consulted for further details. To keep the pa¬ 
per self-contained and before introducing the fermions, 
we will summarize some essential points for the current 
setting. 

In the semi-classical approximation, the Yang-Mills 
partition function is assumed to be dominated by an in¬ 
teracting ensemble of instanton-dyons (anti-dyons). For 
inter-particle distances large compared to their sizes - or 
a very dilute ensemble - both the classical interactions 
and the one-loop effects are Coulomb-like. At distances 
of the order of the particle sizes the one-loop effects are 
encoded in the geometry of the moduli space of the en¬ 
semble. For multi-dyons a plausible moduli space was 
argued starting from the KvBLL caloron [23 that has a 
number of pertinent symmetries, among which permuta¬ 
tion symmetry, overall charge neutrality, and clustering 
to KvBLL. Since the underlying calorons are self-dual, 
the induced metric on the moduli space was shown to be 
hyper-Kahler. 

Specifically and for a fixed holonomy A4(oo)/2u;o = 
i/t 3 /2 with w 0 = nT and r 3 / 2 being the only diago¬ 
nal color algebra generator, the SU(2) KvBLL instanton 
(anti-instanton) is composed of a pair of dyons labeled by 
L, M (anti-dyons by L, M) in the notations of [3]. Gener- 
ically there are N c — 1 M-dyons and only one twisted 
L-dyon type. For the SU(2) gauge group used for most 
of our discussion, M carries (electric-magnetic) charges 
(+,+) and L carries (—,—), with fractional topological 
charges v m = v and vi = 1 — u, respectively. Their corre¬ 
sponding actions are Sl = ‘ZTtVm/ots and Sm = 2,Trvi/a s . 
The M-dyons are also referred to as BPST dyons, while 
the L-dyons as Kaluza-Klein dyons. 

With the above in mind the SU(2) grand-partition 
function is written as 
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Here x m i and y n j are the 3-dimensional coordinate of the 
i-dyon of m-kind and j-anti-dyon of n-kind. Here G[x] a 
(. Kl+Km ) 2 matrix and G[y\ a (Ki+K^) 2 matrix whose 
explicit form are given in mm- V D fi is the streamline 
interaction between D = L, M dyons and D = L,M an- 
tidyons as numerically discussed in [5]. For the SU(2) 
case its Coulomb asymptotic is |T] 


that the dyonic plasma should be dense enough to pro¬ 
duce sufficiently large screening masses, as discussed in 
detail in [Tj. In practice, this limits its application to the 
confined phase with T < T c . The model starts to get 
inapplicable at high density when the dyons are close to 
the maximal packing density. Another limitation is that 
at small enough T the action per dyon 8n 2 /g 2 (T)N c be¬ 
comes small lthereby invalidating the use of the semiclas- 
sical approximation. Our estimates in [J show that the 
model can still be used with reasonable accuracy in the 
range 0.5 T c <T < T c . 


B. Quark zero modes 

For quarks in the fundamental color representation, the 
squared Dirac equation in an external chromo-magnetic 
B and chromo-electric E field, takes the generic form in 
the chiral representation m 


v dd( x - v) -> - 


c D 

a s T 


+ 


I Xm - Vm\ | XL - Vz\ I XM - Vl I I XL - Um I 


( 2 ) 


The strength of the classical Coulomb interaction in |2| 
is Cp/a s = 2.46/a s . At intermediate distances V D p is 
characterized by a core a D p ~ 1/T. The key new ele¬ 
ment in the partition function ([l]) in comparison to our 
previous work [lj , is the introduction of the fermionic de¬ 
terminant detT(x,y) that we will discuss further below. 

The fugacities /,; are related to the overall dyon den¬ 
sity. The dyon density rip could be extracted from lattice 
measurements of the caloron plus anti-caloron densities 
at finite temperateure in unquenched lattice simulations. 
Following [T] we define 


np _ 
T3 “ 



(3) 


with C a constant whose value depends on the regulariza¬ 
tion scheme of the divergent determinant, and ultimately 
on the specific definition of Aqcd- For definiteness, we 
will use 


7T 10, / T \ , A . 

a s 3 ln VO.36 T c ) ^ 

where 10/3 = 11^/6 — AT//3 for N c = 2 and Nf = 
1. The constant inside the logarithm has been fitted to 
lattice measurements of the instanton density for N c = 2 
and Nf = 0. In principle, it should be modified along 
with T c , as the theory changes, e.g. Nf = 0 to Nf = 

1. Since we do not have such lattice data, we will only 

modify the beta function coefficient in front. 

We conclude this section by addressing some limita¬ 
tions of the model described by ([!]). One limitation is 


(-V 2 +4S.(BtE))/ =0 (5) 

with iV = id + A and S a the SU(2) spin generators. 
The signs in ([5]) are commensurate with chirality. In the 
absence of spin, scalar quarks do not admit zero modes 
as —V 2 is semi-positive. With spin, zero modes may 
occur when the spin contribution is negative in ([5]) to 
balance the semi-positive scalar contribution. For a self¬ 
dual object B = E and only the negative chirality quark 
can produce a zero mode state through the “magnetic 
moment term” 


(—V 2 + 4 cr • B) tpp = 0 (6) 

In the dyon the last term is er-B « a-f/p 2 at the core size 
p ss 1/vuq. The first term in © is the squared kinetic 
energy. It is bounded by the uncertainty principle and 
of order 1/p 2 . These two terms can balance by an anti¬ 
aligned spin contribution at the core. Hence the negative 
chirality state is bound in the dyon, while the positive 
chirality state is bound in the anti-dyon. 

The explicit fermionic zero modes of the KvBLL in¬ 
stanton were discussed in m- It was noted that for 
large holonomies the zero mode is localized on one of the 
constituent dyon. The zero modes on individual SU(2) 
dyons were made explicit in [l5]. For SU(2) in the center 
disordered phase with v = 1/2 and zero winding, the M- 
dyon only supports fermionic zero modes with periodic 
boundary condition. In the string gauge for the M-dyon, 
the Dirac equation for the zero mode reads asymptoti¬ 
cally 


(<?4 =F ivojQ — ia ■ d) ip^ = 0 (7) 

with pr 1 the upper and lower components respectively. 
The solutions are hedgehogs with ip ± = ■ 

ff ± (r), 
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(d r + 2/r + vu 0 =F 0//3) / ± (?’) = 0 


( 8 ) 


with x = wo r - Since Tr(</?}<p 2 ) = Tr(f?i 772)1 we may sub¬ 
stitute © into (|10| . The Fourier transform of the result 


is 


Normalizable solutions require <j>//3 < ruling out the 
anti-periodic boundary condition. On the other hand, 
the L-dyon zero mode follows from the M-dyon by switch¬ 
ing v to (1 — v) and then performing a time-dependent 
gauge transformation U = e~ tLU ° X4T3 . Thus the L-dyon 
zero mode is anti-periodic for v = 1/2. 

The L- and M-zero modes associated with the higher 
winding sectors labeled by n on R 3 x S 1 jT2] follow from 
the substitution v —> v + n in ([8]). For (f> = (2m + l)7r, 
there are normalizable M-zero modes if — 2ri + 3 < m < 
2ra 4 ~ 1 . For cj) = 2m'ir, there are normalizable M-zero 
modes if — 2 ”/ l ~ 1 < m' < n . The corresponding L- 
zero modes follow by gauge transformation. Thus anti- 
periodic L-zero modes in a sector n correspond to pe¬ 
riodic M-zero modes in the same sector and vice versa. 
Throughout, only the n = 0 sector will be discussed for 
clarity. 


C. Fermionic determinant 

The main issue discussed in this paper is the behavior 
(pairing or collectivization) of the fermionic zero modes 
into what is called in the literature the “Zero Mode Zone” 
(ZMZ). The approximations used in its description fol¬ 
lows closely the construction, developed for instantons 
and described in detail in refs [14| . The fermionic deter¬ 
minant can be viewed as a sum of closed fermionic loops 
connecting all dyons and antidyons. Each link - or hop¬ 
ping - between L-dyons and L-anti-dyons is described by 
the elements of the “hopping chiral matrix” T 


T(p) = ‘f(\A 1 (p )\ 2 + \A' 0 (p )\ 2 (12) 

with 


A n (p) = 


V 2 w 


. n +h 


dx x 


n+ 


Wn 


i sin(pa;) th| 
px i/shx 


(13) 


with p = p/ujq. The two integrals in (131 for n = 0,1 will 


be carried numerically for our results below. 

Since we will be interested in the center symmetric 
phase of the dyon-anti-dyon ensemble we have left out the 
repulsive linear interaction between unlike dyons (anti¬ 
dyons) in the KvBLL instanton as each sector already 
assume them dissociated. This pair interaction acts as a 
linearly confining force in the center disordered phase. A 
schematic and local pair-interaction between dyons and 
anti-dyons induced by the fermions will be added to © 
below. 

Depending on the dyon density and locations, the de¬ 
terminant can either be dominated by small (binary) 
loops, or very long loops connecting macroscopically 
large number of dyons. The first phase is called “molec¬ 
ular” and is dominated by dyon-anti-dyon clusters, rem¬ 
iniscent of the molecules in the instanton ensemble [231 . 
The second phase of very long loops is called “collec¬ 
tivization” and leads to a nonzero quark condensate. 


D. Bosonic fields 


T(£,2/)-( _ t ^. q J ^ (9) 

with dimensionality (K L + Ki) 2 . Each of the entries in 
T ij is a “hopping amplitude” for a fermion between an 
L-dyon and an L-anti-dyon, defined via the zero mode 
of the dyon and the zero mode (pjj (of opposite chirality) 
of the anti-dyon 


Following mi the moduli determinants in |T]) can 
be fermionized using 4 pairs of ghost fields x'l Xl. m 
for the dyons and 4 pairs of ghost fields miXl,m f° r 
the anti-dyons. The ensuing Coulomb factors from the 
determinants are then bosonized using 4 boson fields 
i>l,Mi w l,m for the dyons and similarly for the anti- 
dyons. The result is a doubling of the 3-dimensional free 
actions obtained in 1 


T ^ = T (x.i - yj) = J d 4 z <p\y(z - Xi)i (7 ■ d)p D (z - %Q 10 ) 

The exact zero modes in the hedgehog gauge relate to the 
zero modes in the string gauge through r/Aa = -~<pf € / 3 a 
with indices A for color and a for spinors. They are 
explicitly given by m 


T]Aol — 


^0 tli 2 
2-\/87r v^a/sha; 




( 11 ) 


Si F [x,v,w] = J d 3 x 

(|Vxl| 2 + I v Xm| 2 + • Vwl + V vm ■ Vwm) + 

(I^Xil 2 + IVxmI 2 + + Vujg • Vwjg)(14) 

For the interaction part V d q, we note that the pair 
Coulomb interaction in ([TJ) between the dyons and anti- 
dyons can also be bosonized using standard tricks [21] [22] 
in terms of a and b fields. We note that a and b are the 
un-Higgsed long range U(l) parts of the original magnetic 
field Fij and electric potential A 4 (modulo the holonomy) 
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respectively. As a result each dyon species acquire addi¬ 
tional fugacity factors such that 

M:e~ b - ia L:e b+ia M : e~ b+ia L : e b ~ ia (15) 


with G a 2 x 2 chiral block matrix 


G = 


0 G(x,y)\ 

-G (x,y) 0 ) 


( 21 ) 


Note that these assignments are consistent with those 
suggested in 0E] using different arguments. As a result 
there is an additional contribution to the free part (14) 


with entries TG = 1. The Grassmanian source contribu¬ 
tions in ( 20 ) generates a string of independent exponents 
for the L-dyons and L-anti-dyons 


S 2 f[ct, 6] = — J d 3 x (V5 • V5 + Ver • Vcr) (16) 

and the interaction part is now 

Si[v,w, b,a,x\ = - J d 3 x 

e -b+icrf M (47 TVm + \ Xm - Xl \ 2 +v M - v L ) e WM ~ WL + 
e +b ~ ia fL (47 tvi + | xl - Xm\ 2 +v l - v m ) e WL ~ WM + 
e~ b ~ la fM ( 4 KVrn + | Xm ~ XL? + v m ~ vz) e Wst ~ WL + 
e + b +ia f L (4^-+ \ XL + V L - Vfj) e WL ~ Wa (17) 


Kl 

gXltV'iGid+blGi.dxi 
i —1 
K L 

1=1 


( 22 ) 


The Grassmanian integration over the X i in each factor 
in ( 22 ) is now readily done to yield 


n \-^\^i{x L i)\ n (23) 

* 3 


without the fermions. We now show the minimal modi¬ 
fications to ( |17| ) when the fermionic determinantal inter¬ 
action is present. 


for the L-dyons and L-anti-dyons. The net effect of the 
additional fermionic determinant in rtd) is to shift the 
L-dyon and L-anti-dyon fugacities in (|17|) through 


E. Fermionic fields 

The determinant for the hopping fermionic zero mode 
can be fermionized using standard methods. For that, 
each entry T (x — y) in (JT|) can be viewed as a cross two- 
body dyon-anti-dyon hopping matrix with a two-body 
inverse TG = 1. To fermionize the determinant we 
define the additional Grassmanians x = {XiiX 2 ) T with 
i,j = 1 ,~,K l L and 


/l -1 = -/i'0 t 7+V> 

ft 


(24) 


where we have now identified the chiralities through -)± = 
(l± 75 )/ 2 . The fugacities f M Xl are left unchanged since 
they do not develop zero modes. 


F. Resolving the constraints 


det T 


= / D[x\ e xfTx 


(18) 


We can re-arrange the exponent in (18) by defining a 
Grassmanian source J(x) = (Ji(x), J 2 (x)) T with 


In terms of ( T4fl7 ) and the substitution (241, the dyon- 
anti-dyon partition function (jTJ) for Nf = 1 can be ex¬ 
actly re-written as an interacting effective field theory in 
3-dimensions, 


k l 


M X ) = '22 Xl$ 3 (x - XLi) 


i =1 
Kf 


M x ) = ^ 22 x 2 & 3 {x-yij) 

3= 1 


(19) 


and by introducing 2 additional fermionic fields ip{x) = 
{ip 1 (x),ip 2 {x)) T . Thus 


e x+T x = I D H’\ exp (- / ^G ip + / JV + J ^ J) ^ 


J dD[ip\ exp (— f f/dG ip) 


Z X [T\ = J D[ip] D[ X ] D[v\ D[w\ D[a] D[b] 

xg-ftp-ftF-Sj-s^ ( 25 ) 

with the additional Nf = 1 chiral fermionic contribution 
= ip^Gip. In the presence of the fermionic fields ip 
and the screening fields <7, b the 3-dimensional effective 
field theory (25) is not integrable. Simple approximation 


schemes will be developed to address this effective action. 


Note that the effective action in (25) is linear in the 
v m l,m,L- These are auxiliary fields that integrate into 
delta-function constraints. However and for convenience, 
it is best to shift away the b, a fields from |l7|) through 
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wm — b + ia^r wm 
w m - b-ia —► wm 


(26) 


which carries unit Jacobian and no anomalies, and re¬ 
cover them in the pertinent arguments of the delta func¬ 
tion constraints as 


-j-V 2 wm + f M e WM ~ WL 
47r 

—/l^ I 7+'0 e WL ~ WM = ^V 2 (b-ia) 

47r 

-j-V 2 w l — f M e WM ~ WL 

47r 

+/ L ^ t 7 +^e“ i -“ M =0 (27) 

and similarly for the anti-dyons. To proceed further the 
formal classical solutions to the constraint equations or 
wm,l [ct, 6 ] should be inserted back into the 3-dimensional 
effective action. As in [5 we observe that the classical 
solutions to (f27|) can be used to integrate the w's in (25) 


to one loop. The resulting bosonic determinant cancels 
against the fermionic determinant after also integrating 
over the y's in (251. The result is 


Zi[T\= J D[i/>]D[<r]D[b]e 

with the 3-dimensional effective action 


S F [a,b] + 

j cPxip^Gil) 

-47r/ M n m 

/A,e— +e - 

+4:TTf L Vl J 

d 3 x^ 1+ ^e WL ~ WM 

+4?r f L vi J 

d 3 x^j-ipe WL ~ w ^ 


(28) 


(29) 


Here S F is S-jf in (161 plus additional contributions re¬ 


sulting from the wm,l (cr, b) solutions to the constraint 
equations ( |27| ) after their insertion back. This procedure 
for the linearized approximation of the constraint was 
discussed in [Tj for the case without fermions. 


III. SU(2) QCD WITH ONE QUARK FLAVOR 


The ground state is parity even so that Jl,m = Jl m 


and 


£5 = 0. By translational invariance, the SU(2) ground 
state corresponds to constant a , 5, £. The classical solu¬ 
tions to the constraint equations (271 are also constant 


(31) 


(e WM ~ WL ) 0 = V/lS/ 2 / m 

and similarly for the anti-dyons. 

A. Effective potential 


The effective potential V for constant fields follows 
from (301 by enforcing the delta-function constraint (30) 


and parity 


- V/V 3 = -MAE 

+^TTf M Vm (e WM ~ WL - 
+2 irf L viZ(e WL ~ WM 


« Wm-Wt 


(32) 


with V 3 the 3-volume. For fixed holonomies v m j, the 
constant w's are real by (271 as all right hand sides van¬ 
ish, and the extrema of (32) occur for 


e w M -w L = ± \ZYif L v l /2f M v r 

e«=±v^/2M 


(33) 


(33) are consistent with (31) only if vi = v m = 1/2 and 
vi = Vm = ~l/2. That is for confining holonomies or a 
center symmetric ground state. Thus 


- V/V 3 = iA£ + 8 tV/l/m£/2 


(34) 


We note that for £ = 0 there are no solutions to the 
extrema equations. The holonomies are no longer con¬ 
strained to the center symmetric state. Since £ = 0 
means a zero chiral condensate (see below), we conclude 
that in this model of the dyon-anti-dyon liquid with light 
quarks, chiral symmetry restoration and the loss of center 
symmetry occur simultaneously. 

For the vacuum solution, the auxiliary field A is also a 
constant. The fermionic fields in (30) can be integrated 
out. The result is a new contribution to the potential 


(34) 


To analyze the ground state and the fermionic fluctu¬ 
ations we bosonize the fermions in (28) by introducing 


two delta-functions and re-exponentiating them 

Z,[T\ = J D[i/j\ D[a} D[b\ D[£] 7J[£ 5 ] D[X] D[ A 5 ] 

x e ~S+f d 3 x iA(i/’ t V'+S)+/ d 3 x i\ 5 (ipUj 5 ifi+i: 5 ) 

(30) 


- V/V 3 -A A£ + 8 ttv// l / m £/2 
d 3 p 


TO 


The saddle point of (35) in £ is solution to 


a a 

iX T —-— — 0 —y A — —— 

\/£ ^£ 


( 36 ) 
























7 


after the substitution A —► iX with a = 47 t-\//l/m/2 . 


Inserting (36) into the effective potential (35) yields 


- V/V 3 = y + f ^ln (1 + A 2 T 2 (p)) (37) 


A ' J (27r) 3 
The saddle point for A is 
d 3 p 


a 

2A 


A 2 T 2 (p) _ 

= ro 


(2?r) 3 1 + A 2 T 2 (p) 


(38) 


It is readily checked that (38) enforces the true minimum 


condition d(V/V 3 ) = 0. From (13) we note that AT(p) « 
Xujq/p 6 falls rapidly with momentum for p > p max with 
p^ ax = WqVA- A simple solution to (38) follows from the 
condition AT(0) 1, i.e. 


Vn 


: p 3 
rmax 


= WoVA 


(39) 


The precise value of Vo is not important as it will be 
traded for the dyon density below. Note that for the 
opposite case of AT(0) <lwe have Vo ~ A 2 /w 0 - This is 
the dilute dyonic density limit which is not our case. The 
dyon ensemble in the center symmetric phase is dense . 
In terms of (391 all equations can be solved analytically. 


However, we have checked that their accuracy is limited. 
All the analysis to follow will be carried exactly without 
these estimates. 


related the integral we called Vo to the dyonic density 
no- So given no, the solution to the gap equation (44) 


fixes A and thus the quark constituent mass M(jp) and 
therefore, through the delta-function constraint in (30), 
the value of E. 

In our approach the np> is not an external input, but 
should itself be calculated from the derivatives of the free 
energy, e.g. the M-dyon density is 


_ 1 S(-V/V 3 ) 

HM 2 d\nf M 


d 3 p M 2 (p) 

(27t) 3 p 2 + M 2 (p) 


= V 0 (44) 


Since np> = Um + til + n M + n L = 4 tim as all partial 
dyonic densities are equal in the confined phase, we have 
Vo = n D / 4. 

TM(p)/A 



p/T 


B. Gap equation 

The free energy depends on two parameters, the mean 
values of A and E fields, which should be chosen at the 
minimum of it. The equations following from vanishing 
first derivatives are known in literature as the “gap equa¬ 
tions”. It is useful, to recast it in terms of the integral 
Vq defined above. In particular, we have 


r _ 4V 0 2 


2Vo 

A 


(40) 


while the effective potential (37) is 


-V/V 3 = 2Hi + /|El„(l + ^) (41) 

We have introduced the momentum-dependent con¬ 
stituent quark mass M(p ) as 


FIG. 1: The momentum dependent quark constituent mass 
TM(p )/A versus p/T. 


C. Chiral condensate 

The non-vanishing of E signals the non-vanishing of 
the chiral condensate (qq) and therefore the spontaneous 
breaking of chiral symmetry . Standard demonstration 
of that is done via introduction of a nonzero but small 


light quark mass to, which changes (23) to 


W-^u) + m] -il’ifoiyLj) + m] (45) 


A rerun of the bosonization scheme with (45) shows that 


only one contribution in (34) is now shifted 


a 


M(p) = ApT(p) = —pT(p) 
ZVq 


(42) 


which is seen to vanish linearly at p/w o ^ 1 and as 1/p 2 
for p/ujq 1. In Fig. [l]we show the behavior of dimen¬ 
sionless mass ratio TM(p )/A as a function of p/T. (42) 
through (38) obeys the gap equation for the A parameter 


d 3 p M 2 (p) 


(27t) 3 p 2 


■ M 2 (p) 


n D 

~T 


(43) 


W/l/mH/2 -»• W/ l / m (E/2 + to) 


(46) 


changing the saddle point solutions (36) to 


A = 


y/S + 2 TO 


and (38) to 


(47) 








































5W mA = 


d 3 p A 2 T 2 (p) 


(48) 


(2tt) 3 l + A 2 T 2 (p) 

The effective potential is now 

- V/V 3 = y + 2mA + J ^ In (l + A 2 T 2 (p)) (49) 


Inserting (491 in the general definition of the chiral con¬ 


densate in the saddle point approximation 


M 

T 


d(V/V 3 ) 

dm 


and using the gap equation we obtain 


M 

T 


= -2A 


(50) 


(51) 


We have used that a is independent of m and that the 
contribution multiplying dX/dm is zero thanks to the gap 
equation. In the chiral limit, A is fixed by the solution of 
the gap-equation ( |43| ) for the constituent quark mass and 
the dyon density tid- It is therefore an implicit function 
of ud , he. A = A [no]. 

The use of (|3j) into (43) leads to only numerical results 
for A and thus M(p). In Fig. [2] we show the behavior 
of the absolute value of the quark condensate | (qq) \/T 3 
versus the dyon density nn/T 3 . Note that j (qq) \/T 3 
decreases with decreasing dyon density. In this range, a 
best fit gives 


m 

T 3 




1.63 


(52) 


For analytical estimates, we note that we can always se¬ 
lect a temperature in the range 0.5 < T 0 < T c for which 
the chiral condensate is \qq\/T ) 3 = 1 at that tempera¬ 


ture. From (51) we have Aq = T 2 /2 in the chiral limit. 


Inserting this value for M(p) in ( |43| shows that the cor¬ 
responding dyon density at T = T 0 is uq/Tq =0.80. 


D. Screened Polyakov lines 


As we noted earlier in (33), the spontaneous break¬ 


ing of chiral symmetry with a finite value of S/2, still 
preserves center symmetry with Vi = v m = 1/2. How¬ 
ever, strict confinement is lost because of screening. 
Heavy fundamental color charges are now screened by 
the light constituent quarks through the formation of 
tightly bound heavy-light and colorless mesons. The 
bound mesons are blind to the Z 2 center and thus to 
the holonomies, with now 


(L(x)) 


0 -p(Y./2+m+0(ot s )) 


(53) 



n D ir 


FIG. 2: The absolute value of the (dimensionless) quark con¬ 
densate | (qq) \/T 3 versus the (dimensionless) dyon density 
n D /T 3 . 


The 0(a s ) contribution in (53) is UV sensitive and re¬ 
quires a specific subtraction. Using (51) together with 
(40) and ((431) where N c = 2, we can recast (53) in the 


chiral limit into the generic relation 


n .d ~ N c (qq) In ((L(x))) (1 + 0(a s )) 


(54) 


for the dyonic density in the range 0.5 < T < T c . (54) 


provides for an independent estimate of the dyon den¬ 
sity in unquenched QCD. Finally, we note that for large 
separation, the correlation of two Polyakov lines clusters 


(L\x)L(0))*\(L(0))\ 2 *e 


—/3(E+2 m+0(a s )) 


(55) 


with a vanishing of the electric string tension due to light 
quark screening. 


IV. MESONIC SPECTRUM 

The stability of the vacuum solution with Nf = 1 can 
be tested by fluctuating in the fermionic channel which 
consists of both a scalar er meson and a pseudo-scalar 
7/ meson. Both are massive. The former through the 
spontaneous breaking of chiral symmetry with finite E, 
while the latter through the U^(l) anomaly with a finite 
topological susceptibility. 


A. Sigma meson 

A simple way to probe the scalar spectrum is to note 
that in the spontaneously broken state, the fermion ki¬ 
netic contribution in (27) is now 


0 G(x,y) 
-G (x,y) 0 


*A1 
-G (x,y) 


xy 


G (x,y) 

^A1 XV 


(56) 
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with 1 xy = S 3 (x — y) .The scalar meson in the long wave¬ 
length limit can be identified with the fluctuations in the 
chiral condensate through iX in (56) or 


zX —y zX 0 T zSX = zXq [IT —t— 

Js 


(57) 


with Ao = — ( qq ) /2 T in the chiral limit. Fluctuations 
in A induce also fluctuations in E. Thus, consistency 
requires 


<5E 


(58) 


Inserting (571 into (56) and (581 into the delta-constraint 

;he effective action for the 


allow for a derivation of t: 
fluctuating parts 5 E, S A. The linear contributions are 
zero by the saddle point equations. So the net contri¬ 
butions are quadratic and higher. In leading quadratic 
order 6E, S A mix, 


A. 




d 3 q (M +9 _TM_ 9 +) 2 
(2tt) 3 (M% + q%)(Ml +q 2 _) 


(64) 


Here, M± = M(q±) and p± = q±p/ 2, where M(q) is 
the running constituent mass in (42). For Nf = 1 we 
have mixing between the scalar as a qq quark state and 
the scalar glueball. The no /2 contribution in (64) is just 
the mixing contribution, while the second contribution 
is clearly the qq quark bubble contribution. A similar 
mixing in the scalar sector was observed in the instanton 
liquid model of the QCD vacuum [22] . 

A comparison of the small momentum expansion of 


(62) after subtraction of the no /2 glue-mix, yields the 
canonical scalar action in x-space 


S 2 [o- s l = 


~2T 


d 3 a 


(|Vcr s | 2 + m 2 s a s 2 ) 


(65) 


S 2 [<5A, AS] = -- 


d 3 p 
(2t r) 3 
d 3 p ri£) 


with 


(2tt)3 4Eq 
d 3 p 


(—2i)5T,(p)SX(—p) (59) 

<SE(p)<JE(-p) 


(2t r) 3 


SX(p)G s (p)SX(-p) 


with 


Gs(P ) = 


d 3 q 2(—Aq T G(< 7 2 )G(p T q) 2 ) 


( 2 tt ) 3 (G 3 (q 3 ) + Xl)(G 3 ( q +pY) + Xl) 


(60) 


To undo the mixing in (60) we solve for <5E to leading 
order 


2i5X - ^<5E = 0 
2 L 0 


(61) 


and insert it back into (60) to give finally the quadratic 


mlft = T 


d 3 q 4 q 2 M 2 {q) 

(27t)3 (q 2 T M 2 (q)) 2 


( 66 ) 


In Fig. [ 3 ] we show the behavior of m 2 / 2 /T 4 from (66) as 
a function of the scaled dyon density no/T 3 . Increasing 
density amounts to lower temperature, with a transition 
density expected around 1. 

mlfhr* 



action for the scalar meson 


n D /r 


FIG. 3: rn 3 f 3 /T 4 versus the dyon density nn/T 3 . 


with 


S2[cr s ]=T^ J -^J^a s (p)A+{p)(7 s (-p) 


A +(P) = AqG s 1 (p) T n D 


(62) 


(63) 


The kernel in (631 can be further reduced using the gap 


equation for no- Thus 


B. eta' meson 

To generate the effective quadratic action for the rj' 
meson, we need to fluctuate asymmetrically around the 
chiral condensate in (561 


iX -T «A+ = ?An 1 ± 


IT] 

VWr, 


(67) 
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A re-run of the preceding arguments for the scalar meson 
yields 


S 2 M = -^2 J V (. p ) A -{ pM ~ p ) ( 68 ) 

with 


A-(P) = ~ 


d 3 q M + M_ (M + M _ + q + ) 
(2+i 3 (gi+M2)(gl+Mi) 


(69) 


for arbitrary current mass to. Using the gap equation for 
njj for non-zero to, we may further reduce (69) into 


n D 


A _(p) = 2 toA + —— + 


1 f d 3 q ( q+M_ - g_M + ) 2 

27 (2^)3 (4 + M 2 )(<£ + Ml ) 


(70) 


Since TA_(0) = / 2 ?n 2 it follows that 


= -m < gg > + Xt 


(71) 


with Xt = TYid/ 4. The first contribution in (75) is the 
Gell-Mann-Oakes-Renner contribution to the q'mass as a 
would-be Goldstone boson, while the second contribution 
is a Witten-Veneziano contribution-like. It suggests an 
unquenched topological susceptibility of \t — Trip/A as 
opposed to the quenched topological susceptibility [J. In 
the chiral limit with m = 0, 


fn m 


Tn 




: XT = 


D 


(72) 


Since the topological susceptibility for dyons Xt ~ 
G(N°) and / 2 « 0{N C ), the rf mass is seen to vanish 
at large N c . In Fig. [4] we display the ratio of the squared 
scalar to pseudo scalar mass as a function of the scaled 
dyon density n d /T 3 assuming f s « f v by chiral symme¬ 
try. The ratio decreases with increasing dyon density or 
lower temperature. This ratio may be compared to the 
value in the QCD vacuum, i.e. T = 0, N c = 3, Nf = 2+1, 
which is about 1/2. Going in the opposite direction, to 
smaller densities, note that the expected phase transition 
density is around 1. 

A simple but approximate understanding of ( f7T| ) fol¬ 
lows by noting that (30) is U(l)y symmetric but upsets 
U(1 )a symmetry through the fermionic contributions. 
Under U(l)^ with ip -+ e' n6S ^ 2 ip the fermionic contri¬ 
butions in (32) change 


+ 2 tt f^e^vi ip^j + ip e WL WM + 27r/^e l6 viip^"f-ip e WL Wm 

(73) 


This amounts to shifting f L i —> f L i 
ity symmetric effective potential (351. 


e ±t8 in the par- 
An estimate of 


m 2 s lm 2 n 



n D IT 3 


FIG. 4: The mass ratio squared, for scalar to pseudo¬ 
scalar mesons m 2 /m 2 , versus the (dimensionless) dyon den¬ 
sity hd/T 3 , assuming f s « f v (see text). 


mass of the if follows by identifying 0/2 -+ r]/f v /\/2 
with a constant rf field. So the rf mass is related to the 
topological susceptibility xt- Specifically, the pertinent 
contribution from the effective potential (35) is now 


- V/V 3 ->• 4tt (/ l / m (E/2 + to)) 1/2 cos (V 2 r)/f v J (74) 

where we have retained the small quark mass m. Ex¬ 
panding (74) yields the quadratic contribution 


T f to 2a 

2 [V^f 2 


2<uv / e'\ 2 

~ir) v 


(75) 


Recall that a = 47T\//m/l/ 2. From ( |5l| ) we have 
2 a/Yi = — (qq) and from (40) we have ayE = 2Vo = 
no/ 2- The squared rf mass follows as in but with 
an incorrect 4yr contribution. 


V. DYON PAIRING THROUGH FERMION 
EXCHANGES 

Above the temperature of the chiral phase transition, 
T > T x , there is no quark condensate, and dyons and 
anti-dyons pair into neutral “dyon-antidyon molecules” 
bound through fermion exchanges m , a situation some¬ 
what reminiscent of the BKT transition dig 

This can be seen by noting that the chiral matrix T 
in Q is banded with a band range set by the inverse 
temperature 


Tjj ~ t f e~^° —>■ t f Sij (76) 

With increasing temperature, the range of T is reduced to 
the nearest neighbor. As a result, the hopping is stalled 
with 


det T 


I k l +k l 


s k l k l 
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The light quark spectrum is now gapped at A± = ±|t/| 
with a vanishing chiral condensate (qq) = 0. An estimate 
of the hopping parameter follows from 


t f = 


d 3 p 

(2tt) 3 


T {p) « 0.8 w 0 


(77) 


using ( TofTT > and Parseval equality. 

A simple but crude estimate of the transition density 
at which the pairing into molecules overtake the chirally 
broken phase is when the molecular gap becomes larger 
than E, thus restoring chiral symmetry. E characterizes 


the size of the delocalized zero mode zone. Using (771 this 
occurs for |t/| ss 2.51T ss E. Since E/T = riDj\\qq) |, 
this means a transition when | (qq) \/nn ~ 1/2.5. From 
the numerical fit (521 this estimate yields to a chiral 
restoration for a diT 


ute dyon ensemble with 


n D 

T 3 


^<^« 0.16 


(78) 


Near the transition temperature, a substantial amount of 
dyons can already be paired, resulting in a weakening of 
the chiral condensate. 

In terms of (76) the partition function 0 is highly 


correlated. The result after summing over pairs is 


Z moi [T\ = j D[b] D[a} D[ X ] D[w} e~ s °~ SM 

x\t f \^F^ L I 2 (\t f \^F~ L ) (79) 

with So defined in ( I4p6 l and Sm defined in fll7| for 
only M and M. The argument of the modified Bessel 
function I 2 is composed of 


F l = 


F- l = 


j d 3 xe +b ~ ia+WL - WM 
x f L (Attvi + \xl - Xm | 2 + V L - v M ) 

J d 3 x e +b+ia+WL ~ wgI 

x /l ( 4 fvj + |XL — Xm 1 “ + V L ~ v m) 


(80) 


The molecular partition function in (79) is highly non¬ 
linear in the auxiliary fields. For large |tj|, we may use 
the asymptotic form of I 2 (z) « e 2 /v / 27r z in (80) and 
linearize the argument of the modified Bessel function 


y/Ftfi 


Fl + F- l 

V2 


1 - 


n 


F i 


2 (F l + F- l Y 


(81) 


As a first step to be justified below, we may drop the 
non-linear contributions in (81) and the pre-exponent in 


(79) to have 


Z mo i[T]tt J D[b]D[a]D[x\D[w] e ~ s oSM+(FL+F L )/V 2 


(82) 


after rescaling f hL -)> f L ,l\ t f\/V2 in F. (82) is 
now analogous to the SU(2) Yang-Mills partitionlunc- 
tion [DEI with the new re-scaled fugacities. A rerun of 
the arguments in this case shows that the ground state 
is still center symmetric. However, the ground state is 
chirally symmetric. It is worth noting that this state 
is parity even, so the neglected non-linear corrections in 
(81) amounts to (1 —1/4) which is about a 25% reduction 
in the pertinent pressure contribution which is then 


F mol = 


lnZ mo \ //m/l |t/|\ ^ / co n 

WF~ SnT {-^-) (83) 


VI. HIGHER NUMBER OF COLORS AND 
FLAVORS 


The extension of the current analysis to many N c col¬ 
ors and Nf massless flavors is straightforward in princi¬ 
ple. For finite N c the KvBLL instanton splits into N c 
constituent dyon with 1/N C topological charge and fu- 
gacity fi with 1 < l < N c . The L-dyon zero mode which 
is anti-periodic is now carried by the l = N c constituent 
dyon. The net effect is a change in the fermionic contri¬ 
bution in (30) through G —> G ® 1/, and a change in the 


parity even effective potential (32) as 


- V/V 3 -a +i\N f E + AnfiPi ( e Wi ^~ Wi + e w '+^~ w ') 

+4?rf Nc v Nc -^det Nf (V’b+V’s) e WNc ~ WNc+1 

+4 FfN c v N c ^det Nf (^7-^s) e w ^~ w ^+i 

(84) 


where the implicit i-summation is over i = 1,..., N c — 1. 
The potential V has manifest SU (Nf)y x SU(A r /)^ x 
U(l) v flavor symmetry. As a result, the parity even 
effective potential (35) after pertinent bosonization and 
Fierzing yields 


-V/V 3 —> +i\NfT, + 2a(N c ) E® (85) 

+n, !wp ln (i-XT 2 M) 

with x = Nf/N c and a(N c ) = Anf(N c )/2 x . The mean 
fugacity is 


f(N c ) = (fi....f Nc ) 1/Nc (86) 
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We note that its scaling with N c follows from the scaling 
of each fugacity by semi-classics, i.e. /) ss 1/a 2 . Thus 


f(N c ) « N 2N ^ w IV 2 


so that a(N c ) « TV 2 . 


(87) 


d 3 p 

(27t) 3 


< 90 > 


Using the exact form (13) and the solution to the gap 


equation at T = T 0 , we have 


A. Gap equation and chiral condensate 

For general x = Nt/N c , the saddle point equation in 
S of ((851) gives 


E = 


A 


2xa(N c ) 


( 88 ) 


after the shift —iA —> A and A = Nf A. With this in mind 
and inserting (881 into (85) yields 


d 3 p rr \'2 r ^ 10-37 

7y%T \p) = —- 

(27 r) 3 To 


(91) 


which shows that (901 is in general upset, and this case 


does not possess a minimum. 

With this in mind and for x < 2, the extremum of (89) 
in A yields the gap-like equation 


A 


\2 xa(N c ) J 

with the new identification 


NcVnj- i 
a(N c ) 


(92) 


- V/V 3 = -2 a(N c ) (x - 1) 


+xN c 


2 xa(N c ) 


(89) 




V Nf -i = 


d 3 p 


N- 


^-(p) 


(2n) 3 i + ^ T 2(p) 


(93) 


The dyonic density is now identified with 


The case x = 1 is special. The effective potential in (85) 
is linear in E with no a priori saddle point along E. We 
have checked that taking the saddle point in A first, and 
then the saddle-point in E after the substitution results in 
the same gap equation to follow. Also, it can be checked 
explicitly that the same results follow by taking the limit 
x —> 1. 


The effective potential (891 has different shapes de¬ 
pending on the ratio of the number of flavors to the 
number of colors x. Let us explain that in details for 
four cases: 

(i) If x < 1 the first term in (89) has a positive coef¬ 
ficient and a negative power, so it is decreasing at small 
A. At large value of A the second term is growing as InA. 
Thus a minimum in between must exist. This minimum 
is the physical solution we are after. 

(ii) If 1 < x < 2 the coefficient of the first term is 
negative but its power is now positive. So again there is 
a decrease at small A and thus a minimum. 

(iii) If x > 2 the leading behavior at small A is now 
dominated by the second term which goes as A 2 with 
positive coefficient. One may check that the potential is 
monotonously increasing for any A with no extremum. 
There is no gap equation, which means chiral symmetry 
cannot be broken in the mean-field approximation. 

(iv) If x = 2 there are two different contributions of 
opposite sign to order A 2 at small A. An extremum forms 
only if the following condition is met 


no = = 2A r c^ / -i (94) 

2 olnjM 


In terms of ( fOApM | the running constituent mass M ( p ) = 
(\/Nf)pT(p) obeys the gap equation 


d 3 p M 2 (p ) 


(2-7t) 3 p 2 


■ M 2 (p) 


- w < 95 > 


From (88) it follows that A « N c Nf. From (92) it follows 
that Vjv7_i ~ N c and therefore no ~ N 2 . The dyonic 
description we have reached is consistent with large N c 
counting. Since nn ~ N 2 1 crystallization in the form 
of dyonic salt is expected at large N c j26] . 

For x < 2 the center symmetric vacuum also breaks 
spontaneously chiral symmetry, with a vacuum conden¬ 
sate given by 


M 

T 


— — 2An 


(96) 


Ao is the value of A in the chiral limit. Since Ao ~ NfN c 
the chiral condensate in (96) is of order NfN c as ex¬ 
pected. 

To summarize: Chiral restoration as well as the loss 
of center symmetry occur simultaneously for x x > 2 as 
per our result in (34). This value of x = Nf/N c is close 
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to the the critical value of x x = 5/3 originally suggested 
in the instanton liquid model |14| (first reference). First 
simulations of the dyon ensemble with N c = 2 m also 
indicate that the border line seems to be Nf = 4, in 
agreement with x x ss 2. 

Current lattice data are summarized e.g. in Fig.5 of 
[25] for N c = 3. Indeed, they seem to indicate a change 
in the value of the chiral transition temperature (in units 
of the vacuum string tension) T x /y/a at x = 2 or Nf = 
6, but instead of vanishing, this ratio remains flat up 
to N f = 8. For such a large Nf the number of quark 
lines 2Nf connected to an L dyon is large. Maybe the 
correlations between them are too strong for the mean 
field approximation to remain valid. 

The chiral transition we have discussed in this section 
should not be confused with another phase transition in 
theories with a large number of flavors, namely the con¬ 
formal (fixed infrared coupling) phase. The reported lat¬ 
tice m and holographic (Veneziano limit) [25] results 
put this conformal transition at a much larger number of 
flavors ^conformal ~ 4, or Nf = 12 for N c = 3. 


(981 as a function of the reduced dyon density fin at 
intermediate densities. The linearity of the logarithm in 
fiD both at small and asymptotic dyon densities follow 
from the scaling of Vo = nn/2N c with A as discussed in 

©. 

In terms of (98) the classical pressure contribution in 


(971 simplifies to 


Ptot ~ Pp 
N r T A 


(1 - x(l - k)) n D 


( 100 ) 


with 1 < k < 2. In the quenched limit or x = 0 it reduces 
to the dyonic result obtained for the pure Yang-Mills 
analysis in [Ij ignoring the Debye-Huckel corrections. 
The fermion induced interactions in the center symmet¬ 
ric phase increase the pressure away from the free limit 
for 0 < x < 1. Remarkably, for small densities(l — re) « 0 
and (100) with fermions is close to a free ensemble of 


dyons. For large densities or large N c , (1 — k) « —1 and 
(100) with fermions is more repulsive than the free dyon 


ensemble. 


B. Thermodynamic of dyonic phase with x < 1 


In the presence of light quarks, the total thermody¬ 
namical pressure of the dyon-anti-dyon liquid consists 
of the classical and non-perturbative contributions in 
(85) at the extremum, plus its perturbative correction 


for finite and symmetric holonomies v = 1/N C [13] . plus 
the purely perturbative black-body contribution (ignor¬ 
ing the higher order 0(a s ) quantum corrections). Identi¬ 
fying the classical pressure with —V/fiV^ with /? = 1/T, 
we have 


Ptot ~ Ppei 
N C T 4 


+(1 - x) h D 


(97) 



[ d \ln 


! (2tt)3 

V P 2 ) 


with hp = nr>/(N c T 3 ). The assessment of the logarith¬ 
mic integral follows by numerical integration using the 
explicit form of M{p) and the solution to the gap equa¬ 
tion. The result is linear in the reduced dyon density for 
small and asymptotic densities 



n D IT 


3 


FIG. 5: The fermionic loop ln[F]/T 3 versus the reduced dyon 
density tid/T 3 . See text. 


The perturbative contribution is given by 



45 

7tt 2 x 
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1) 


7n 2 x 

180 


N 2 


( 101 ) 


1 

T 3 


d 3 p 
(2t r) 3 


In 



M 2 (p) \ 
P 2 ) 


n(n D )n D 


(98) 


with < 1) « 1 and n{hu 1) ~ 2. A simple 

interpolation to the overall numerical results is 


I -l- 

K ^°) ~ 1 , H (") 

1 ' 10 

Since no ~ 0(N C ), large density corresponds to large N c 
with k~ 2, modulo crystalization. In Fig. [5] we display 


The first contribution is the free gluon contribution in 
the symmetric phase with v = 1 /N c . The second contri¬ 
bution is the free black-body gluon contribution, which 
is cancelled by the second contribution in leading order 
in l/iV c in the symmetric phase [3|. The third contri¬ 
bution is the free quark contribution in the symmetric 
phase with v = 1/N C . The fourth and last contribution 
is the black-body quark contribution, which we note is 
cancelled by the third contribution in leading order in 
1 /N c . This generalizes the observation in [3] to QCD. 

An estimate of the transition temperature T c from the 
symmetric phase with v = 1 /N c to the asymmetric phase 
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with v = 0 follows when all the non-black-body contri¬ 
butions in the total pressure Vtot cancel out. This occurs 
when the rescaled dyon density n D = nn/(N%T 3 ) solves 


no c 



j_\ i + 7 -j 

N%) 1 + x(n{n Dc ) - 1 ) 


( 102 ) 


with again 1 < /«(n£> c ) < 1. 


C. Thermodynamic of molecular phase 

For completeness we note that near the chiral transi¬ 
tion most of the LL dyons start to pair into molecules, 
for which case the total pressure is more appropriately 
described by 


tot, mol 


-v 0 


2A 4 


rp4 


at 


N f 


It/. 

V2 


(103) 


with |ty | = |t/|/A. Here, the scale parameter A = A/T is 
identified with the vacuum dyon de nsity no -> 2A ^ /a 2 s T 
as in [3]. We note that for Nf —> 0, (103) is off by 2“ J in 
the ground state pressure from the Yang-Mills limit in [3]. 
This can be traced back to our linearized approximation 
in (81). The critical temperature is now 


T r 


with a, = N r ~ 


2A 


\i f \ Nf 


V2 J hi{x) 


(104) 


h(x) = 


45 


lx 

T 


l - 


m 


(105) 


(1041 characterizes the transition from a center symmet¬ 


ric but chirally symmetric phase to a center asymmetric 
and chirally symmetric phase. Which likely transition is 
to occur first can be estimated by comparing the total 
liquid pressure in (97) to the total molecular pressure in 


(103). This is best addressed using mixtures. 


VII. CONCLUSIONS 

We have extended the mean field treatment of the 
SU(2) dyon-anti-dyon liquid in [I], to account for light 
quarks. Anti-periodic fundamental quarks develop zero 
modes for the L, L-dyons only. In the dense phase under 
consideration with T < T c , these zero modes are collec¬ 
tivized into a Zero Mode Zone of quasi-zero modes which 
dominates the low-eigenvalue part of the Dirac spectrum. 


This phenomenon is analogous to the one used in the in- 
stanton liquid model }14j , although the zero modes them¬ 
selves and most of the results are different. The impor¬ 
tant interplay between center symmetry and the sponta¬ 
neous breaking of chiral symmetry which is absent in M 
is now clarified. In particular, we have explicitly shown 
how the chiral effective Lagrangian for light quarks knows 
about confinement. 

In the infrared, the fermionic determinant is entirely 
saturated by these quasi-zero modes [T5], modifying the 
dyon-anti-dyon measure initially suggested in Ei to 
include light quarks. For the Nf = 1 case of one mass¬ 
less quark, we have shown that the fermionic determi¬ 
nant modifies the L,L-dyonic fugacities through chiral 
fermionic bilinears that upset the Ua( 1) symmetry. By 
a series of bosonic and fermionic techniques we have ex¬ 
plicitly mapped the interacting dyon-anti-dyon Coulomb 
liquid with light quarks on a 3-dimensional effective the¬ 
ory with fermions. The translationally and parity invari¬ 
ant ground state was shown to follow from pertinent gap 
equations. The ground state breaks spontaneously chi¬ 
ral symmetry by developing a fermion condensate. For 
Nf = 1 it does not produce a Goldstone mode because of 
the 11^(1) anomaly. We have derived explicit expressions 
and estimates for masses of the a and 77 mesons. 

We have shown how the model generalizes to arbitrary 
number of flavors and colors. In the whole temperature 
interval in which our approach is applicable, the ensemble 
is center symmetric (confining) and breaks spontaneously 
chiral symmetry provided x = Nf/N c < x x « 2. The 
loss of center symmetry and chiral symmetry restoration 
in this model seem to occur simultaneously for x > 2. 
We have noted that in the case of a very large Nf the 
fermion-induced interactions maybe too strong to trust 
the mean field approximation we used. This point needs 
to be pursued numerically on the lattice. 

This conclusion can be compared to the critical value 
of x x ss 5/3 of the numerical simulation of the instanton 
liquid model )14l (first reference). The first simulation of 
the dyon ensemble with fermions, for N c = 2 m , found 
the border line case is Nf = 4, also in agreement with 
x x « 2. So whether the transition is an artifact of the 
mean field approximation or not remains to be studied. 

The chiral transition should not be confused with the 
transition to conformal - fixed infrared coupling- phase, 
for which current lattice and holographic results put this 
transition at a much larger number of flavors ^conformal ~ 
4 , or Nf = 12 for N c = 3. 

Near and above the chiral transition the fermionic cor¬ 
relations are strong enough to pair L-dyons with L-anti- 
dyons into molecules. The dilute regime involved has 
been explored numerically in m- In this paper we only 
produced some estimates of the transition parameters. 

In our approach the confining and chirally broken 
phase have been treated via the mean field approximation 
only, so the resulting gap equation has either a finite or 
zero E, with a finite jump. In the future one can prob¬ 
ably include the LL correlations in the ensemble. The 
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result would be a depletion of the chiral condensate with 
perhaps a more continuous cross-over transition as cur¬ 
rently observed in QCD-like theories with several flavors 
of massive quarks. 

In the extreme case where all L-dyons and L-anti-dyons 
pair to molecules, we have shown that the linearized 
molecular partition function supports a phase with cen¬ 
ter symmetry but restored chiral symmetry. It would be 
interesting in the future to see if such a phase may exist, 
at some N c ,Nf. At this moment, lattice data on that 
issue are also not clear, see e.g. [27] , 

An estimate of the pressure in the center symmetric 
phase shows that both the free gluon and fermion loop 
nearly cancel out in leading order in N c . We have used 
it to estimate the transition density from a center sym¬ 
metric phase to a phase with broken center symmetry. A 
similar estimate of the transition density was made in |T] 
in the absence of fermions. 

The current model can be expanded and improved in a 
number of ways. The current analysis has been done for 
the sector with zero 9. As indicated earlier, the L-zero 
modes were selected over the M-zero modes, creating a 
topological unbalance and a lack of manifest 0-periodicity 
in the induced effective action. This can be generalized 
to arbitrary 6 angle through an extended formalism. 


Also, we have not included here some Coulomb correc¬ 
tions discussed in [T], to keep the analysis simpler and 
to illustrate the interdependence of the center symme¬ 
try and the chiral symmetry breakings in this model. 
These corrections are Debye-like and still within the 
semi-classical analysis. 

Some improvement of the moduli space metric may 
be considered in the future. We recall that the moduli 
space metric used in (|T|) while exact for LM dyons at 
all separations, is only exact asymptotically for LL, MM 
dyons. While the formers attract, the latters repel. In 
the center symmetric or confining phase, we expect the 
like-dyons to stay away from each other while the un¬ 
like dyons to mingle and screen. In the dyon-antidyon 
channels the treatment is so far classical only, with one- 
loop effects absent. We hope to report on some of these 
issues, as well as on a full analysis of the meson spectrum 
for the dyon-anti-dyon liquid with Nf > 1 next. 
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